Morphology and scaling in the noisy Burgers equation 
Soliton approach to the strong coupling fixed point 
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The morphology and scaling properties of the noisy Burgers equation in one dimension are treated 
by means of a nonlinear soliton approach based on the Martin-Siggia-Rose technique. In a canonical 
formulation the strong coupling fixed point is accessed by means of a principle of least action in 
the asymptotic nonperturbative weak noise limit. The strong coupling scaling behaviour and the 
growth morphology are described by a gas of nonlinear soliton modes with a gapless dispersion law 
and a superposed gas of linear diffusive modes with a gap. The dynamic exponent is determined by 
the gapless soliton dispersion law, whereas the roughness exponent and a heuristic expression for 
the scaling function are given by the form factor in a spectral representation of the interface slope 
correlation function. The scaling function has the form of a Levy flight distribution. 
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There is a current interest in the scaling properties and 
general morphology of nonequlibrium models. Here the 
noisy Burgers equation in one dimension provides maybe 
the simplest continuum description of an open driven 
nonlinear system This equation has the form of 

a conserved Langevin equation 

du/dt = tA7 2 u + XuVu + Vrj , (1) 

where v is a damping constant characterizing the linear 
diffusive term, A a coupling strength for the nonlinear 
mode coupling term, and r\ a Gaussian white noise driv- 
ing the system into a stationary state and correlated ac- 
cording to 

(r](xt)r](x't')) = AS{x - x')S(t - t') . (2) 

The equation (jl]) has been discussed extensively in partic- 
ular as a model for the stochastic dynamics of the slope 
field u = V/i for the self-affine growth of an interface 
subject to annealed noise from the drive or environment 
described by the Kardar-Parisi-Zhang (KPZ) equation 

dh/dt = vV 2 h+(\/2)\{Vh) 2 +r, . (3) 

The main emphasis has been on the scaling properties 
of the slope field in the large x - long t limit as embodied 
in the dynamical scaling form j|pj 

(u{xt)u{x't')) = \x- x'r 2(1 ~ c) /(|i - t'\/\x - x'\ z ) , (4) 

where £ is the roughness exponent, z the dynamic expo- 
nent, and / the scaling function; f(w) cx u;"^ 1 ^)/ 2 for 



large values of its argument. In this context the nonlinear 
Galilean symmetry of (|l]) 

X — > X — Xugt , U — > U -\- Uq (5) 

with A entering as a structural constant implies the scal- 
ing law 

( + z = 2. (6) 

Another property specific to one dimension is the sta- 
tionary probability distribution [|],(|||] > 

P(u) cx exp [-(v/A) J dxu 2 } . (7) 

independent of A. Hence u is an independent Gaussian 
variable, (u(x)u(x')) — (A/2v)5{x — x'), and the height 
field h = J dxu performs random walk with exponent 
( = 1/2. The scaling law (||) then yields z — 3/2, defin- 
ing the Burgers-KPZ universality class. This result also 
follows from dynamic renormalization group calculations 
HQ but the associated infrared stable strong coupling 
fixed point does not seem accessible by the methods of 
renormalized perturbation theory and the e-expansion. 

In the linear Edwards- Wilkinson (EW) case for A = 
the slope correlations are 

Ak 2 

{u{ku)u(-k - u)) = ^ 2 ~ {vk2)2 , (8) 

yielding the exponents £ = 1/2 and z = 2, characterizing 
the EW universality class with scaling function 
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f(w) = (A/2^)(47r^)~ 1/2 w^ 1/2 exp [-1/Avw] 



(9) 



In a recent letter [0 we approached the strong cou- 
pling fixed point behaviour by means of a mapping of 
the Burgers equation (P onto a solid-on-solid model and 
further onto a discrete spin 1/2 chain model. In a har- 
monic oscillator representation valid for large spin com- 
bined with a continuum limit this approach eventually 
yields a Hamiltonian description and a set of coupled 
equations of motion with soliton solutions with dispersion 
E oc p 3 / 2 , yielding the exponent z — 3/2, characterizing 
the strong coupling fixed point. The approach further- 
more gave a tentative picture of a growing interface in 
terms of a soliton gas of paired solitons representing a 
growing step in the height field. 

Here we present a unified approach to the noisy Burg- 
ers equation based on a Martin-Siggia-Rose path integral. 
This method provides a generalization of the stationary 
distribution (jTj) to the time-dependent nonlinear case and 
yields the following physical picture and scaling results: 

1) In the weak noise limit the general morphology of a 
growing interface in the stationary regime is described by 
a Landau-type quasi-particle many-body theory in terms 
of a gas of nonlinear soliton modes with superposed lin- 
ear diffusive modes. 

2) The noise-driven fluctuations correspond to "quan- 
tum fluctuations" in the underlying non-Hcrmitian re- 
laxational "quantum description" , yielding the quasi- 
particle picture. 

3) The scaling exponents and the scaling function follow 
as a by-product from the soliton dispersion law and the 
spectral representation of the correlations. The many- 
body formulation explains the robustness of the rough- 
ness exponent under a change of universality class. 

From a field theoretical point of view we identify the 
noise strength A in (||) as the effective small parameter. 
Furthermore, the fundamental probability distribution 
or path integral has an essential singularity for A = 0. 
Hence our approach is based on a nonperturbative sad- 
dle point or steepest descent approximation to the path 
integral. It is precisely in this respect that the dynamic 
renormalization group method based on an expansion in 
A and in the noise contraction A fails to access the strong 
coupling fixed point. 

The special role of the noise is easy to understand. The 
turning on of the noise in the Burgers equation is a sin- 
gular process in the sense that even weak noise A ~ 
will eventually drive the system into a stationary state 
whereas for A = the noiseless deterministic Burgers 
equation (assuming vanishing currents at the boundaries) 
describes a damped interface, although with pattern for- 
mation or dissipative structures formed by the nonlinear 
cascade term in the transient time regime. In the pres- 
ence of noise the system passes from the transient to the 
stationary regime at a crossover time t co determined by 
the noise strength. In the linear EW case t co is of or- 



der l/(^fc 2 ) log (1/A) showing that t co — > oo for A — > 0; 
k is the wave number of the particular mode considered 
(note that t co — > oo also for k — > 0, characteristic of a 
conserved hydrodynamical mode). In the linear EW case 
the singular character of the noise strength A is reflected 
in the stationary distribution (]?]) which has an essential 
singularity for A = 0. 

It is well-known that the noiseless Burgers equation for 
A = HHH, du/dt = v\I 2 u + AuVu, supports parity- 
breaking right hand solitons connected by ramp solutions, 
corresponding to cusps connected by convex parabolic 
segments in the height field. Superposed on the solitons 
are linear diffusive modes with a gap in the spectrum. 
This configuration constitutes the growth pattern in the 
transient time regime. For A = the solitons disappear 
and the equation reduces to the linear diffusion equations 
supporting gapless diffusive modes. 

In the stationary noise-driven case for A ^ the soli- 
ton excitations and diffusive modes also play an impor- 
tant role in the analysis of the physics and scaling be- 
haviour of the strong coupling fixed. In the noisy case 
the solitons appear as stationary or saddle point solu- 
tions to the path integral in the weak noise limit; the dif- 
fusive modes correspond to Gaussian fluctuations about 
the stationary points. Also the noise fluctuations lift the 
broken parity symmetry and excites both right and left 
hand solitons. As regards the diffusive modes both de- 
caying and growing modes are excited in accordance with 
the time reversal invariance of the stationary regime as 
for example reflected in the evenness in uj of the slope 
correlations (||) in the linear case. 

In terms of the Martin-Siggia-Rose path integral Q| 
the slope correlations are given by Bj 



{u{xt)u{x't')) = 



I Ylxt dudcp exp [i{v/ A) S]u{xt)u(x' t 1 ) 



JUxt dud<pexp[i(v/A)S] 



(10) 



which provides the generalization of the stationary EW 
case (u(x)u(x')) = jW x duP(u)u(x)u(x r ) with P(u) 
given by (^). The action S in (|o|) has the form 



S 



dxdt[udip I dt — Ti.(u, ip)] 



(11) 



with Hamiltonian density 

H = -i{v/2)[{\/uf + (V<yj) 2 ] + (A/2)w 2 V<^ . (12) 



The correspondence of (|l0l,|llj,|12[) with the Feynman 
phase space path integral || allows us to consider Ajv 
as an effective Planck constant and the slope field u and 
noise field if as canonically conjugate variables satisfy- 
ing the Poisson bracket {u(x), <p(x')} — 5(x — x'). It is 
also characteristic of a stochastic growth problem that 
the Hamiltonian, being the generator in the relaxational 
master equation, in general is complex. The imaginary 
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harmonic part drives the diffusive modes whereas the real 
nonlinear term describes the intrinsic growth of the inter- 
face. The path integral replacing the noisy Burgers equa- 
tion is deterministic, the noise-driven fluctuations being 
replaced by the contributing paths or configurations in 
the path integral. 

The canonical form of the path integral allows for a 
principle of least action. In other words, in the weak 
noise limit A — > we can determine "classical" solu- 
tions or orbits from the stationary points of the ac- 
tion. The weak noise solutions are then given by the 
equations of motion, du/dt — {H,u} = —SH/Scp and 
dip/dt = {H,tp} = Sn/5u, 



du . 2 



\ Y7 ^ 

XuVu , — 
at 



+ivV 2 u - 



XuVip . 
(13) 



The deterministic nonlinear coupled field equations Jl3| ) 
for the slope and noise fields thus replace the noisy Burg- 
ers equation (Q) in the asymptotic nonperturbative weak 
noise limit. We note that i) the equations are invariant 
under a shift of <p reflecting the conserved noise in (^) and 
ii) like the noisy Burgers equation, invariant under the 
Galilean transformation (^) supplemented with <p — ► cp. 

In addition to the constant solutions (u, (p) = (uq, ipo), 
corresponding to the zero-energy stationary state, the 
field equations also support nonlinear soliton solutions of 
both parities. The static solitons are given by 



±u+ tanh (k s (x — xq)) 



k s = Xu+/2v . (14) 



Here k s is a wavenumber depending on the soliton am- 
plitude, setting the inverse length scale, and xq a center 
of mass position. Note that only the right hand soliton 
appears in the noiseless case. Boosting the soliton to a 
finite velocity v and denoting the boundary values by u± 
the general soliton condition is || 



u + + w_ 



-2v/\ 



(15) 



relating the soliton amplitude and off-set to the propa- 
gation velocity. 

For A = the field equations yield diffusive modes 
u = A exp (—iujt + ikx) + c.c. with gapless dispersion 



= —ivk 



(16) 



The Hamiltonian (|l2) is harmonic, the path integral Jic| ) 
Gaussian, and it is an easy task to evaluate the slope cor- 
relations (H). 

It is a well-known feature of field theoretical saddle 
point soliton or instanton calculations that multi-soliton 
solutions also contribute to the stationary points ||. 
Here this implies that the extremal path in the weak noise 
limit corresponds to a dilute soliton gas of right and left 
hand solitons matched according to the soliton condition 



(jl5|). In other words, the path integral formulation di- 
rectly produces a physical picture of the morphology of a 
growing interface. The natural variable is the slope field 
u. A single soliton configuration connecting two station- 
ary slope configuration corresponds to the propagating 
top or bottom of a growing step in the height field h. A 
single moving step in h is represented by a pair of solitons 
matched according to Jl5| ) . In a similar fashion a moving 
plateau in h corresponds to four solitons, a growing tip to 
a three-soliton configuration, etc. The morphology of a 
growing interface can thus be interpreted as a dilute gas 
of solitons corresponding to moving steps in the height 
field. We note that the noise radically changes the mor- 
phology from the noiseless transient case in that the noise 
excites both "up" and "down" cusps in h. Superposed on 
the soliton gas are linear diffusive "ripple" modes with a 
gap in the spectrum. 

The canonical formulation of the path integral enables 
us to associate energy, momentum, and action with a 
given soliton configuration. This allows for a dynami- 
cal selection criterion similar to the Boltzmann factor 
exp (— E/T) in equilibrium statistical mechanics which 
associates an energy E with a given configuration con- 
tributing to the partition function; in the dynamical case 
the action S provides the weight function for the dynam- 
ical configuration. 

The energy is given by E = J dxTi whereas the 
momentum follows from the Poisson bracket algebra, 
P = J dxuVtp. For a single soliton we have in terms 
of the boundary values, E = ±i(\/6)(u 3 i _ — v£_) and 
P = ±i(l/2)(u?_ — u 2 J). In the case of a single moving 
soliton satisfying w_ = we obtain in particular 



£=(4i/3)M 3 /A 2 



P 



-(2i)v\v\/X 2 



(17) 



We note the nonlinear velocity dependence characteristic 
of soliton solutions [ |io|]ll) . Eliminating v we obtain the 
soliton dispersion law E = A(\/2/3) exp (i9)P 3 / 2 , where 
the phase 6 = tt + (7r/4)(u/|u|). 

The linear diffusive mode spectrum in the presence of 
the solitons is analyzed by a linear stability analysis of 



the field equations (13). Like in the noiseless case |6 the 
associated eigenvalue problem is exactly soluble [12 . In 
addition to zero-eigenvalue bound states corresponding 
to the soliton translation modes lifting the broken trans- 
lational symmetry, the spectrum also has a band of phase 
shifted diffusive modes with a gap in the dispersion law 



lo = -iv(k 2 + k 2 ) 



(18) 



depending on the soliton amplitude according to (|l4|) . 

In order to complete the analysis of the noisy Burg- 
ers equation it remains to discuss the contribution of 
the Gaussian fluctuations in the path integral about the 
solitons. The basic information lies in the path integral 
and will be discussed elsewhere. Here we take a heuris- 
tic point of view and discuss the fluctuations in terms 
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of the underlying non-Hermitian relaxational "quantum 
field theory" leading to the path integral by the usual 
Feynman procedure The canonical fields u and ip 
must now be interpreted as "quantum operators" satisfy- 
ing a commutator algebra, the Hamiltonian ( |12|) is a non- 
Hermitian operator, and the field equations (|f3|) are the 
associated Heisenberg equations of motion. The "classi- 
cal" localized soliton solution (Jf4|) becomes a delocalized 
"quantum" quasi-particle with frequency Q = E/(A/v) 
and wavenumber K — P/(A/v). For a pair of matched 
quasi-particles representing a growing step the "quan- 
tum" dispersion law is 



(u(x,t)u(0,Q)) cx 



iAt, 



(19) 



Evaluating the group velocity of a wave packet repre- 
senting the "quantum soliton" we find the classical prop- 
agation velocity, in complete accordance with an effec- 
tive correspondence principle. In a similar way the linear 
diffusive modes are "quantised" becoming quasi-particles 
corresponding to "ripple" modes on the solitons with dis- 
persion law @ flO). 

The picture of the growing interface is thus complete. 
The noise-induced fluctuations admits a "quantum me- 
chanical" interpretation. The dominant morphology con- 
sists of a fluctuating "quantum" gas of solitons with a 
gapless dispersion law with exponent 3/2. Superposed on 
the soliton gas are linear diffusive "ripples" modes with a 
gap in the dispersion depending on the soliton amplitude. 
In the limit of vanishing growth for A — > the solitons 
disappear and the diffusive modes become gapless cor- 
responding to thermodynamic fluctuations governed by 
the EW equation with Boltzmann weight (Q). 

In the "quantum interpretation" the slope correla- 
tion function (u(xt)u(00)) is given by a time-ordered 
Green's function (0\Tu(xt)u(00)\0) f§. Here u is the 
time-dependent "quantum" slope field and |0) denotes 
the appropriate stationary zero-energy state. Assuming 
t > 0, using the Hamiltonian and momentum operators 
to "displace" the field at (x,t) to (0,0), and inserting 
a complete set of momentum "quasi-particle" states, we 
derive an effective spectral representation which allows us 
to discuss the scaling properties in simple terms, 



u(xt)u(00)) = / dKG(K)e~ i( - nt+Kx) 



(20) 



Here G(K) is an effective form factor which in general is 
complex owing to the non-Hermitian Hamiltonian (12). 
At long distances (20) samples the small wave number 
region in the integral and assuming that the form factor 
is regular for small K and a general dispersion law of the 
form 



= A + Rf- 



(21) 



characterized by the gap A and the exponent [A, we ob- 
tain, absorbing i in a rescaling of K 



x- 1 J dKe -iK»it/*»)-iK . (22) 



First, in the presence of a gap there is no scaling due 
to the exponential prefactor in (p2|), consequently, only 
gapless excitations contribute to the scaling behaviour. 
For A = comparing with (0) we infer z = fx and 
C - 1/2. 

In the linear EW case the diffusive gapless modes 
with dispersion ( |l6| ) exhaust the spectrum and we ob- 
tain fi = z = 2, corresponding to the EW universality 
class. The spectral form (^2|) yields the scaling function 
in ©. 

In the nonlinear Burgers-KPZ case the nonlinear soli- 
ton modes with gapless dispersion ( |l9|) exhaust the bot- 
tom of the spectrum and yields fj, = z = 3/2 for the 
Burgers-KPZ universality class. The linear diffusive 
modes develop a gap, become sub-dominant and do not 
contribute to the scaling. This mechanism explains the 
change in universality class: The universality class is de- 
termined by the dominant gapless excitation. The spec- 
tral form (]2(]) also elucidates the robustness of the ex- 
ponent £ which is the same for both universality classes. 
For the stationary fluctuations we set t — in (20) and 



the resulting scaling function yielding £ does not depend 
on the specific quasi-particle dispersion law. 

Finally, we also obtain a heuristic expression for the 
scaling function 



dKe 



-i(K z w+K) 



(23) 



which incidentally has the same form as the probability 
distribution for Levy flights [jl3| . 

In this letter we have presented the outline of a novel 
approach to the growth morphology and scaling be- 
haviour of the noisy Burgers equation in one dimension; 
details will appear elsewhere. Using the Martin-Siggia- 
Rose technique in a canonical form we have demonstrated 
that the so far elusive strong coupling fixed point be- 
haviour is associated with an essential singularity in the 
noise strength and can be accessed by appropriate field 
theoretical soliton techniques. 

Discussions with J. Krug, M. Kosterlitz and A. Svane 
are gratefully acknowledged. 
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